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Convergence of an Euler scheme for a hybrid stochastic-local 
volatility model with stochastic rates in foreign exchange markets 
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Abstract 

We study the Heston-Cox-Ingersoll-Ross++ stochastic-local volatility model in the context of 
foreign exchange markets and propose a Monte Carlo simulation scheme which combines the 
full truncation Euler scheme for the stochastic volatility component and the stochastic domestic 
and foreign short interest rates with the log-Euler scheme for the exchange rate. We establish 
the exponential integrability of full truncation Euler approximations for the Cox-Ingersoll~Ross 
process and find a lower bound on the explosion time of these exponential moments. Under a 
full correlation structure and a realistic set of assumptions on the so-called leverage function, we 
prove the strong convergence of the exchange rate approximations and deduce the convergence of 
Monte Carlo estimators for a number of vanilla and path-dependent options. Then, we perform 
a series of numerical experiments for an autocallable barrier dual currency note. 
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1 Introduction 

The class of stochastic-local volatility (SLV) models have become very popular in the financial 
sector in recent years. They contain a stochastic volatility component as well as a local volatility 
component (called the leverage function) and combine advantages of the two. According to Ren et 
al. [35], Tian et al. m and van der Stoep et al. [38|, the general SLV model allows for a better 
calibration to European options and improves the pricing and risk-management performance when 
compared to pure local volatility (LV) or pure stochastic volatility (SV) models. We focus on the 
Heston SLV model because the Cox-Ingersoll-Ross (CIR) process for the variance is widely used in 
the industry due to its desirable properties, such as mean-reversion and non-negativity, and since 
semi-analytic formulae are available for calls and puts under Heston’s model and can help calibrate 
the parameters easily. The local volatility component allows a perfect calibration to the market 
prices of vanilla options. At the same time, the stochastic volatility component already provides 
built-in smiles and skews which give a rough calibration, so that a relatively flat leverage function 
suffices for a perfect calibration. 

In order to improve the pricing and hedging of foreign exchange (FX) options, we introduce 
stochastic domestic and foreign short interest rates into the model. Assuming constant interest 
rates is appealing due to its simplicity and does not lead to a serious mispricing of options with 
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short maturities. However, empirical results m have confirmed that the constant interest rate 
assumption is inappropriate for long-dated FX products, and the effect of interest rate volatility 
can be as relevant as that of the FX rate volatility for longer maturities. There has been a great 
deal of research carried out in the area of option pricing with stochastic volatility and interest rates 
in the past couple of years. For instance. Van Haastrecht et al. m extended the model of Schobel 
and Zhu [36| to cnrrency derivatives by including stochastic interest rates, a model that benefits 
from analytical tractability even in a full correlation setting due to the processes being Gaussian, 
whereas Ahlip and Rutkowski [T], Grzelak and Oosterlee [18] and Van Haastrecht and Pelsser im 
examined the Heston-2CIR/Vasicek hybrid models and concluded that a full correlation structure 
gives rise to a non-affine model even under a partial correlation of the driving Brownian motions. 

There has also been an increasing interest in the calibration to vanilla options of models with 
stochastic and local volatility and stochastic interest rate dynamics. For instance, Guyon and 
Labordere |20| examined Monte Carlo-based calibration methods for a 3-factor SLV equity model 
with a stochastic domestic rate and discrete dividends, while Cozma, Mariapragassam and Reisinger 
m and Deelstra and Rayee |13) examined 4-factor hybrid SLV models with stochastic domestic 
and foreign rates and proposed different approaches to calibrate the leverage function. Hambly, 
Mariapragassam and Reisinger |2I] took a different path and discussed the calibration of a 2-factor 
SLV model to barriers. As an aside, in this paper, the spot FX rate is defined as the number of 
units of domestic cnrrency per unit of foreign currency. 

The model of Cox et al. [TO] is very popular when modeling short rates because the square- 
root (CIR) process admits a unique strong solution, is mean-reverting and analytically tractable. 
Cox et al. found the conditional distribution to be noncentral chi-squared and Broadie and Kaya 
[8] proposed an efficient exact simulation scheme for the square-root process based on acceptance- 
rejection sampling. However, their algorithm presents a number of disadvantages such as complexity 
and lack of speed, and it is not fit to price strongly path-dependent options that require the value 
of the FX rate at a large number of time points. Furthermore, in the context of a stochastic-local 
volatility model, the correlations between the underlying processes make it difficnlt to simulate a 
noncentral chi-squared increment together with a correlated increment for the FX rate and the 
short rates, if applicable. 

As of late, the non-negativity of the GIR process is considered to be less desirable when modeling 
short rates. On the one hand, the central banks have signihcantly reduced the interest rates since 
the 2008 financial crisis and it is now commonly accepted that interest rates need not be positive. 
On the other hand, if interest rates dropped too far below zero, then large amounts of money would 
be withdrawn from banks and government bonds, putting a severe squeeze on deposits. Hence, we 
model the domestic and foreign short rates nsing the shifted GIR (CIR-I—1-) process of Brigo and 
Mercurio [7]. The GIR-t-1- model allows the short rates to become negative and can fit any observed 
term structure exactly while preserving the analytical tractability of the original model for bonds, 
caps, swaptions and other basic interest rate products. 

In the present work, we put forward the Heston-Gox-Ingersoll-Ross-I—I- stochastic-local volatil¬ 
ity (Heston-2GIR-t-l- SLV) model to price FX options. Independent of the correlation structure, 
the model is non-affine and hence a closed-form solution to the European option valuation problem 
is not available. Finite-difference methods are popular in finance and when the evolution of the 
exchange rate is governed by a complex system of stochastic differential equations (SDEs), it all 
comes down to solving a higher-dimensional partial differential equation (PDE). This can prove 
to be difficnlt due to the curse of dimensionality, because the number of grid points required in¬ 
creases exponentially with the number of dimensions. Monte Garlo algorithms are often preferred 
due to their ability to handle path-dependent features easily and there are numerous discretization 
schemes available, like the simple Euler-Maruyama scheme, see, e.g., Glasserman m However, 
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there are several disadvantages of this discretization, such as the fact that the approximation pro¬ 
cess can become negative with non-zero probability. In practice, one can set the process equal to 
zero when it turns negative - called an absorption hx - or reflect it in the origin - referred to as a 
reflection fix. An overview of the Euler schemes considered thus far in the literature, including the 
full truncation scheme, can be found in Lord et al. [32]. 

The usual theorems in Kloeden and Platen m on the convergence of numerical simulations 
require the drift and diffusion coefficients to be globally Lipschitz and satisfy a linear growth 
condition, whereas Higham et al. |23| extended the analysis to a simple Euler scheme for a locally 
Lipschitz SDE. The standard convergence theory does not apply to the CIR process since the square- 
root is not locally Lipschitz around zero. Consequently, alternative approaches have been employed 
to prove the weak or strong convergence of various discretizations for the square-root process. 
Deelstra and Delbaen |12j . Alfonsi |2], Higham and Mao |22| and Lord et al. |32j examined the 
strong global approximation and found either a logarithmic convergence rate or none at all. Strong 
convergence of order 1/2 of the symmetrized and the backward (drift-implicit) Euler-Maruyama 
(BEM) schemes was established in Berkaoui et al. |6| and Dereich et al. |14j . respectively, albeit 
in a very restricted parameter regime for the symmetrized scheme. Alfonsi [3| and Neuenkirch 
and Szpruch |33| recently showed that the BEM scheme for the SDE obtained through a Lamperti 
transformation is strongly convergent with rate one in the case of an inaccessible boundary point, 
while Hutzenthaler et al. m established a positive strong order of convergence in the case of an 
accessible boundary point. 

Hutzenthaler et al. |28| identified a class of stopped increment-tamed Euler approximations for 
nonlinear systems of SDEs with locally Lipschitz drift and diffusion coefficients and proved that 
they preserve the exponential integrability of the exact solution under some mild assumptions, 
unlike the explicit, the linear-implicit or some tamed Euler schemes, which rarely do. However, the 
results of Hutzenthaler et al. do not apply to the present work beacuse the diffusion coefficient in 
the square-root model is not locally Lipschitz. In this work, we first prove that an explicit Euler 
scheme for the CIR process retains exponential integrability properties, which plays a key role 
in establishing the boundedness of moments of Euler approximations for SDE systems with CIR 
dynamics in one or more dimensions. 

To the best of our knowledge, the convergence of Monte Carlo algorithms in a stochastic-local 
volatility context has yet to be established. Higham and Mao [22] considered an Euler simulation 
of the Heston model with a reflection fix in the diffusion coefficient to avoid negative values. They 
studied convergence properties of the stopped approximation process and used the boundedness of 
payoffs to prove weak convergence for a European put and an up-and-out call option. However, 
Higham and Mao mentioned that the arguments cannot be extended to cope with unbounded pay¬ 
offs. We work under a different Euler scheme and overcome this problem by proving the uniform 
boundedness of moments of the true solution and its approximation, and then the strong conver¬ 
gence of the latter. The existence of moment bounds for Euler approximations of the Heston model 
is important for deriving strong convergence |53| and it has been a long-standing open problem 
until now. Furthermore, the existence of moment bounds for hybrid Heston-type stochastic-local 
volatility models plays an important role in the calibration routine HH. 

In this paper, we focus on the Heston stochastic-local volatility model with CIR-I—I- short interest 
rates and examine convergence properties of the Monte Carlo algorithm with the full truncation 
Euler (FTE) discretization for the squared volatility and the two short rates, and the log-Euler 
discretization for the exchange rate. We prefer the full truncation scheme introduced by Lord et 
al. because it preserves the positivity of the original process, is easy to implement and is found 
empirically to produce the smallest bias of all explicit Euler schemes and to outperform the quasi- 
second order schemes of Kahl and Jackel |29| and Ninomiya and Victoir |34) . We also choose the 
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FTE scheme over the BEM scheme since the latter works only in a restricted parameter regime that 
is often unrealistic, since the BEM scheme is equivalent to the implicit Milstein method employed 
by Kahl and Jackel to first order in the time step, and since the quanto correction term in the 
dynamics of the foreign short rate would lead to technical challenges in the convergence analysis. 

The quadratic exponential (QE) scheme of Andersen [1] uses moment-matching techniques to 
approximate the noncentral chi-squared distributed square-root process and is an efficient alterna¬ 
tive to the FTE scheme that typically produces a smaller bias at the expense of a more complex 
implementation. Empirical results suggest that the QE scheme outperforms the FTE scheme in 
the case of a severely violated Feller condition [39], whereas when the Feller condition is satisfied, 
the two schemes perform comparably well because the bias with the FTE scheme is small even 
when only a few number of time steps are used per year. Since there is no natural extension to the 
QE scheme to multi-dimensional Heston-type models that involve several correlated square-root 
processes such as the ones studied here, we prefer the ETE scheme in the subsequent convergence 
analysis. 

In summary, to the best of our knowledge, this paper is the first to establish: (1) the exponential 
integrability of the ETE scheme for the CIR process; (2) the boundedness of moments of order 
greater than 1 of approximation schemes for Heston-type models; (3) the strong convergence of 
approximation schemes for models with stochastic-local volatility (with a bounded and Lipschitz 
leverage function) and stochastic CIR++ short rates; (4) the weak convergence for options with 
unbounded payoffs, in particular, for European, Asian and barrier contracts (up to a critical time). 

The remainder of this paper is structured as follows. In Section]^ we introduce the model and 
discuss the postulated assumptions and an efficient calibration. In Section we first define the 
simulation scheme and discuss the main theorem. Then, we investigate the uniform exponential 
integrability of the full truncation scheme for the square-root process and prove convergence of the 
exchange rate approximations. We conclude the section by establishing the convergence of Monte 
Carlo simulations for computing the expected discounted payoffs of European, Asian and barrier 
options. In Section we carry out numerical experiments to justify our choice of model and to 
demonstrate convergence. Finally, Section contains a short discussion. 


2 Preliminaries 

2.1 Model definition 

In its most general form, we have in mind a model in an FX market, for the spot FX rate S, the 
squared volatility of the FX rate v, the domestic short interest rate and the foreign short interest 
rate . Unless otherwise stated, in this paper, the subscripts and superscripts “d” and “/” indicate 
domestic and foreign, respectively. Consider a filtered probability space (H,T', P) and suppose that 
the dynamics of the underlying processes are governed by the following system of SDEs under the 
domestic risk-neutral measure Q: 

'dSt = {rf - r{)Stdt + a{t, St)VmStdWi, Sq > 0, 
dvt = k{6 - vt)dt + ^y/vt dW^, vq > 0, 
rt = 9t + hd{t), 

' dgf = Udd - gf)dt + id\[7Uwt, do > 0, 

4 = gl + hf{t), 

Jg{ = {kfdf - kfg( - PsfCfcr{t, St)^/vtg{)dt + ^f^J~g(dW/, > 0, 
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where are standard Brownian motions, a is the leverage function and and 

hf are deterministic functions of time. The mean-reversion parameters k, kd and kf, the long-term 
mean parameters 6, 9d and Of, and the volatility parameters ^d and are positive real numbers. 
The quanto correction term in the drift of the foreign short rate in (2.1) comes from changing 
from the foreign to the domestic risk-neutral measure [9]. As an aside, for Hull-White short rate 
processes, changing from the domestic spot measure to the domestic T-forward measure leads to 
a dimension reduction of the problem because the diffusion coefficients of the short rates do not 
depend on the level of the rates m- Since this is not the case for CIR short rate processes, we 
prefer to work under the domestic spot measure. Note that the above system can collapse to the 
Heston-2CIR-l—I- model if we set a = 1, or to a local volatility model with stochastic short rates if 
we set /c = ^ = 0. The standard Heston SLV model is the special case kd = Cd = kj = = hdj = 0. 

We can also think of (2.1) as a model in an equity market with stock price process S, stochastic 
interest rate and stochastic dividend yield , in which case the quanto correction term vanishes. 
We consider a full correlation structure between the Brownian drivers {W^, W^, W^, W-^}, i.e., no 
assumptions on the constant correlation matrix S are made, where 


1 

Psv 

Psd 

Psf 

Psv 

1 

Pvd 

Pvf 

Psd 

Pvd 

1 

Pdf 

-Psf 

Pvf 

Pdf 

1 


Furthermore, we work under the following assumptions: 


( 2 . 2 ) 


(.4,1) The leverage function is bounded, i.e., there exists a non-negative constant cTmax such that, 
for all t G [0, T] and x G [0, oo), we have 

0 < a{t,x) < amax ■ (2.3) 


(.42) There exist non-negative constants A, B and a positive real number a such that, for all 
t, M G [0, T] and G [0,oo), we have 

\cr{t, x) — a{u, y)\ < A\t — ul^ + B \x — y\ . (2-4) 

(.43) There exists a non-negative constant hmax such that, for all t G [0, T] and i G {d, /}, we have 

\hi{t)\ < hmax ■ (2.5) 


Hence, we assume that a is bounded. Holder continuous in t and Lipschitz in St- According 
to for the leverage function to be consistent with call and put prices, it has to be given by 
the formula (2.9), which depends on the calibrated Dupire local volatility. In practice, the local 
volatility function usually arises as the interpolation of discrete values obtained from a discretized 
version of Dupire’s formula. Hence, there is no loss of generality from a practical point of view in 
assuming that the leverage function is Lipschitz continuous and bounded on a compact subset of 
of the form [0,r] x [xmin,Xmax], and furthermore that 


a(t, x) — a(t A T, Xmin ^X<Xmin 4" ® ^xe{Xmin,Xmax) ^max ^X>Xmax) ' 


( 2 . 6 ) 


Then a is globally Lipschitz continuous and the second assumption holds with a = 1. We also 
assume that hdj are bounded. According to [7], for a perfect fit to the initial term structure of 
interest rates, each of the two shift functions must be given by the difference between the (flat- 
forward) market instantaneous forward rate and a continuous function of time. Then hdj are 
piecewise continuous and the third assumption holds. 
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2.2 Model calibration 


The calibration of the 4-factor Heston“2CIR++ SLV model (2.1) is of paramount importance and 


represents a mandatory step for the efficient pricing of derivative contracts. In HU , we consider the 
4-factor SLV model and propose a new calibration approach that builds on the particle method of 
|20j . combined with a novel and efficient variance reduction technique that takes advantage of PDE 
calibration to increase its stability and accuracy. The numerical experiments in HU suggest that 
this method almost recovers the calibration speed from the associated 2-factor SLV model with 
deterministic rates. We assume here a partial correlation structure where only psd and psf may 
be non-zero, and denote by and the domestic and foreign discount factors associated with 
their respective money market accounts, i.e., for any t G [0,T], 


Df = e~ -^0 


and d( = 6-^0 


rl du 


(2.7) 


In addition to (2.1), we consider a pure local volatility (LV) model, 


dSY^ = {ft - fl)St^dt + ^7Lv(t, St^)St^dWf, = So, 


( 2 . 8 ) 


where for i G {d, f}, ft = — log P^{0, t) is the market instantaneous forward rate at time 0 for a 
maturity t and P*(0,t) is the market zero-coupon bond price at time 0 for a maturity t. 


Suppose that the LV model (2.8) has been calibrated and that ulv has been determined. The 


main building block of the calibration routine is expressing the leverage function a in terms of the 
local volatility function ctlv- We state the necessary and sufficient condition for a perfect calibration 
to vanilla options in the following theorem, proved in HU- 


Theorem 2.1 (Theorem 1 in HU)- If the spot process marginal density function under model (2.1 


is continuous in space and under assumptions (.Al) and {A3), the call price under model (2.1 
matches the market quote for any strike K and maturity T <T* if and only if 


a\T,K) = 


^[Dt^\ST = K] 
¥.[DttvT\ST = K] 


I^) + 


iL2 


a" Cry 

dK‘2 


E 


IT 


ir/,-4){ST-Ky 


- KE 


Dj’{rj’ — ft) I^St>k 


+ KE 


Dt{rt - ft) Ist>k 


(2.9) 


where Clv is the call price under model (2.8), all expectations are under the domestic risk-neutral 
measure, ip = 2 \/2, ( = ^amax and T* is as given below. 


(1) When k < ip(, 


(2) When k > ip(. 


rji^ _ 


[2 


vr 

—h arctan 




T* = oo. 


( 2 . 10 ) 


( 2 . 11 ) 


This result generalizes the formula in m to a stochastic foreign short rate. Deelstra and Rayee 
HU obtained a similar formula for a 4-factor SLV model with Hull-White short rate processes. An 
important step in the derivation is that the stochastic integral 


I^St>K cr{t, Syy/WtDfStdWf 
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is a true martingale. On the one hand, T* is a lower bound on the explosion time of the second 
moment of the discounted spot process DfSt (see Proposition 3.12). On the other hand, moments 
of the Heston model can explode in finite time [5], a property that is inherited by the Heston- 
type model (2.1). Note also Theorems 3.1 and 3.6 for possible moment explosions of the numerical 
approximations. Therefore, the formula (2.9) may not hold for some values of the model parameters 
and for large maturities T. However, in practice, T* is very large. For instance, following our 
calibration routine to EURUSD market data from nn, we found T* = 28.6. 

As an aside, note that if the volatility in model (2.1) is purely local, i.e., if u = 1, we recover 
the formula in [9], whereas if the two short rates are deterministic, i.e., if = 0, we recover 

the formula in |20l [35l, i.e., 


a{T,K) = 


^MT,K) 


( 2 . 12 ) 


I St = K\ 

In the computations in Section we assume a simple correlation structure where only psv may 
be non-zero. We use a five-step calibration routine which is explained in detail in HU. First, we 
calibrate the LV model (2.8) and find ulv- We also calibrate the associated 2-factor SLV model with 
deterministic rates via PDE methods, using formula (2.12) and results from the literature [3135]. 
Independently, we calibrate the two CIR-I—I- processes under their respective markets. Then, we find 
the Heston parameters by calibrating the associated 4-factor SV model, i.e., the Heston-2 CIR-I—I- 
SV model, using the closed-form solution for the call price in pQ, extended to CIR-I—I- processes, 
with other parameters from the above calibration. Einally, we calibrate a in the Heston-2 CIR-I—I- 
SLV model (2.1) by working with the particle method of m and extending it using the calibrated 
2-factor SLV model with deterministic rates as a control variate. For a partial correlation structure 
with possibly non-zero psd and psf, we may use the historical values of the correlations between 
the spot FX rate and the short-term zero-coupon bonds, and extend the approximation formulae 
of [T8j for the call price from a 3-factor Heston-CIR model to a 4-factor Heston-2CIR+-l- model. 


3 Convergence analysis 

3.1 The simulation scheme 

We employ the full truncation Euler (FTE) scheme from [32] to discretize the variance and the two 
short rate processes. Consider the CIR process 

dyt = ky{9y - yt)dt + ^yy/^tdW'l'. (3.1) 

Let T be the maturity of the option under consideration and create an evenly spaced grid 

T = N5t, tn = n5t, Vn G {0,1,..., V}. 

Eirst, we introduce the discrete-time auxiliary process 

yt„+, = yt„ + ky{9y - )(5t + (3.2) 

where y~^ = max(0,y) and and its continuous-time interpolation 


m = yt„ + ky{9y - )(t - tn) + ^yJyt„ - ^1) 


(3.3) 
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for any t £ [tn,tn+i), as suggested in [22]. Then, we define the non-negative processes 


yt = Vt 

yt = vt 


(3.4) 

(3.5) 


whenever t £ [tn,tn+i)- Let v and be the FTE discretizations of v and respectively. Taking 
into account the presence of the quanto correction term in the drift of the foreign short rate, we 
similarly define 


9t=9L+[kf0f-kf{g{y-Psf^fa{tn,St„)^vl{g{yY^-tn)+^f^{g{y{w/-W/j, (3.6) 

where S is the continuous-time approximation of S defined below, as well as 

<al = is!)* (3.7) 

\ s{ = (gL)* (3.8) 

whenever t £ [tn,tn+i)- For i £ {d,f}, the domestic and foreign short rate discretizations are 

rl = 9l + hi{t). (3.9) 

Finally, we use an Euler-Maruyama scheme to discretize the log-exchange rate. Let x and x be the 
actual and the approximated log-processes, and let S = be the continuous-time approximation 
of S. Moreover, define h = hd — hf. Then the discrete method reads; 


ftn + l / ^ f 1 r, - \ 

Xt^+i =xt„+ / h{u)du+ [pt-gl - {tn, St„)vt„jSt + a{tn, 

J tfi 

However, we find it convenient to work with the continuous-time approximation 

xt = xt„+ f h{u)du+(jif - g{ -tn) + 

J tyi 

where AVF/ = Wf — and a{t^St) = a(tn,St„) whenever t £ [tn,tn+i)- Hence, 
xt = XQ + j (j=^-rl-^d‘^{u,Su)vy^du + j d{u, Su)V^dW^. 


(3.10) 


(3.11) 


(3.12) 


Note that the convergence of the continuous-time approximation ensures that the discrete method 
approximates the true solution accurately at the gridpoints. Using Ito’s formula, we obtain 


St = So+ f {r^-rQSudu+ f a{u, Su)\/v^ SudW^. 
Jo Jo 


(3.13) 


We prefer the log-Euler scheme to the standard Euler scheme to discretize the exchange rate process 
because the former preserves positivity and produces no discretization bias in the S'-direction when 
a is constant. Furthermore, if v, g^ and gJ are also constant, then the log-Euler scheme is exact. 






3.2 The main theorem 

Define the arbitrage-free price of an option and its approximation under (3.13), 


t/ = E 




[/ = E 


- f 
3 Jo 




dt 


fis) 


(3.14) 

(3.15) 


The payoff function / may depend on the entire path of the underlying process and the expectation 
is under the domestic risk-neutral measure. 


Theorem 3.1. If2kf6j- > and under assumptions (^1) to (.4.3), the following statements hold. 

(1) The approximations to the values of the European put, the up-and-out barrier call and any 
barrier put option defined in (3.15) converge as 5t —>■ 0. 

(2) ifc = i^. 


and 




f>2fc ) 


(3.16) 


then the approximations to the values of the European call, Asian options, the down-and-in/out 
and the up-and-in barrier call option defined in (3.15) converge as 6t —)• 0. 


Remark 3.2. If assumption (.41) holds, then for the purpose of this paper we choose the smallest 
upper bound on the leverage function, namely 


CTmax = sup {a{t,x) \ t G [0,r], X G [0,oo)}. 


(3.17) 


Remark 3.3. If the exchange rate and the foreign short rate dynamics are independent of each 
other, i.e., if psf = 0, then the quanto correction term in the drift of the foreign short rate vanishes 
and Theorem 3.1 holds independent of the Feller condition 2kf9f > albeit with a simpler proof. 
If the domestic and foreign short rates are constant throughout the lifetime of the option and if. 


moreover, a{t,x) = 1 for all t G [0,T] and x G [0, oo), then the system of equations (2.1) collapses 
to the Heston model and Theorem |3 .1 1 applies with C = ?• This extends the convergence results of 
Higham and Mao [22] to options with unbounded payoffs. 


Remark 3.4. Tian et al. m calibrated the term-structure Heston SLV model parameters to the 
market implied volatility data on EURUSD using 10 term structures and maturities up to 5 years. 
The data in Table suggest that the condition in (3.16) is typically satisfied in practice, both for 
shorter and longer maturities. 


Table 1: The calibrated Heston SLV parameters for EURUSD market data from August 23, 2012 |T7| . 


T 

k 


^max 

C 


1 month 

0.885 

0.342 

1.600 

0.547 

11.8 years 

5 years 

0.978 

0.499 

1.300 

0.649 

9.3 years 


In HH, the Heston-2CIR-l—I- SLV model is calibrated to the EURUSD market data from March 
18, 2016, for maturities up to 5 years. The data in Table indicate that the condition in (3.16) is 
satisfied even for very long maturities. Furthermore, the following parameter values for the foreign 
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Table 2: The calibrated Heston-2CIR++ SLV parameters for EURUSD market data from March 18, 2016 

[in- 


T 

k 


max 

c 


5 years 

1.412 

0.299 

1.399 

0.418 

32.3 years 


short rate process are recovered: kf = 0.011, 6f = 1.166, = 0.037, so that the Feller condition 

also holds, i.e., 2kf9f = 0.0257 > 0.0014 = 0. 

In equity markets, the mean-reversion speed is usually several times greater than the volatility 
of volatility. For instance, Hum et al. [25| calibrated the Heston model for the S&P 500 index from 
January 1990 to December 2011 using a combination of two out-of-the-money options, and found 
that k = 1.977 ;§> 0.456 = Furthermore, we typically have C < Ij so that the condition in (3.16) 
holds even for longer maturities. 


3.3 The square-root process 


In order to prove the convergence of the approximation scheme in (3.13), we first need to examine 
the stability of the moments of order higher than 1 of the actual and the discretized processes. 
However, this problem is directly related to the exponential integrability of the CIR process and its 
approximation. Let y be the CIR process in (3.1) and y be the piecewise constant FTE interpolant 
as per (3.5). 


Proposition 3.5. Let A > 0 and define the stochastic process 

Qt ^ exp J yu du^, Vt > 0. 

IfT< T*, then the first moment of Qt is bounded, i.e., 

E [©t] < oo, 

where T* is as given below. 

(1) When ky < 


rji^ _ 


V2AC2 _ ^2 [2 


TT 

—h arctan 


\/2A^2 _ f,2j_ 


(2) When ky > 


T* = oo. 

Proof. Follows directly from Proposition 3.1 in [5]. 

Theorem 3.6. Let A > 0 and define the stochastic process 

Qt^explxf yudu\, Vt > 0. 


10 


IfT<T* and 6 t < ky^; then the first moment of Qt is uniformly bounded, i.e., 

sup E [Qt] < OO, 


(3.18) 


(3.19) 


(3.20) 

(3.21) 
□ 

(3.22) 

(3.23) 


where T* is as given below. 


10 















(1) When ky < \/0.5A^y, 


(2) When ky > VO.SA^y, 


rjt^ _ 


'\/ 2 A ky 


rp* _ “^ky 


(3.24) 


(3.25) 


Proof. First, we prove that there exists rj > 1 independent of 6t such that, for all uj G [0,1], 

- 2r]UjkyT - 2r/ + 2 < 0. (3.26) 

Fix any r] > 1 and define the polynomial 


frjiui) =u T] XfyT - 2ujr]kyT - 2(r/ - 1), 


with two distinct real roots 


UJ± = 


ky ± kl + 2{r] - l)Xil 

riKlT 


Since a;_ < 0 < w+, we have that f'qiiO, 1]) < 0 if and only if frj{l) < 0, i.e., 

v‘^X(lT^ - 2r?(l + kyT) + 2 < 0. 


(3.27) 

(3.28) 

(3.29) 


This holds for some rj > 1 if and only if the quadratic polynomial in rj on the left-hand side has a 
real root greater or equal to one. Hence, we find the necessary and sufficient conditions: 

{V^Cy-ky)T <1, and 2Xf‘^T < ky + yjkl + AXil ky + < 2XilT < iky, 


which are equivalent to T <T*, with T* defined in (3.24) ~ (3.25). Fix any rj satisfying (3.26). 


Next, we prove by induction on 0 < m < that, for all 6t < 5 t, 

N—m—l 

X5t 


E [Qt] < exp ^{kyOy + Uyf,y) X{6t)'^{m - l)m E 


exp 


N—m— 

\xst y: 


yt, + r]mX6tytj^_^^ 


i=0 


where 


£2 


1 

+ 




^ 47r(l — ky^TY 27r Y 4(1 — kyST^ 




, (3.30) 


(3.31) 


Let , 0 < t < r} be the natural filtration generated by the Brownian motion W^ and consider 
the shorthand notation IE* [ • ] = E [ • |^|^] for the conditional expectation. Note that (3.30) clearly 
holds when m G {0,1}. Next, let us assume that (3.30) holds for 1 < m < A7 and prove the inductive 
step. Conditioning on the cr-algebra we obtain 


E [©r] < exp ^ -v{ky0y + I'yf.y) X{6tY (m — l)m 


X E 


exp 


TV—m—1 

X6t yp 

i=0 


tjV-m-l 


exp ymXStytj^_^^ 


(3.32) 
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For convenience, define w = If ^ ~ AA(0,1), then JL =^/MZ. Let X 

be the inner expectation in (3.32), then 


X < Eq,, 


exp 


^ym\5t max 0, tc + ky{9y — w)6t + ^yVw5tZ | 


There are two possible outcomes, namely re = 0, in which case 

X < exp ^r]mky6yX{6t)‘^'^, 

and w > 0, which is treated now: 


X < 


'^0 




exp 


+ 


/ 


^0 1 


—oo a/ 27r 


-- 2 ;^ + r]m^yX\/w z + rimX6t[w + ky{6y — w)5t] dz 


1 


exp < —-Z > dz 


where 


zo = 


kyOydt + (1 — ky5t)w 

^yVwdt 


Recall that 6t < ky^ and define 


zi = zq — r]m^yXy/w{6t)^^‘^. 

Then zi < zq < 0. If cj) and <h are the standard normal PDF and CDF, then 
I<^( ' ' . 1 

and hence 

where 


a = 2r/m(l — ky6t) + rf'rin?‘^^X{6t)‘^ > 0. 

Applying the mean value theorem to <h G , we can find z G [zi, zq] such that 

^>( 2 : 0 ) - ^{zi) = ( 2:0 - zi) cl) [z] < [zQ - zi) cj) ( 2 : 0 ). 
Hence, from (3.35) and ( 3.36| ), 

<^{zo)-Hzi) < ^r?me,A(5t)3/2.V^exp|-J^^Mi±^^M!H 




2Qw6t 


(3.33) 


(3.34) 

(3.35) 

(3.36) 


( 2 : 0 ) + exp yqmX6t[ky6ydt + (1 — ky6t)w] + -r]^m‘^^yX‘^w{6t)^j^l — <1> (2;i)| 
X < exp i^rimky9yX{6t)^ + -aAr(;(5t||l + 4> ( 2 : 0 ) — 4> ( 2 ;i)|, 


(3.37) 

(3.38) 


(3.39) 


= b, constant w.r.t. w 

We can think of the right-hand side as a function of w, say / : (0, 00 ) 1 —>■ M. Next, we show that 

\2 w„., ^ n (3 40) 


f{w) < r]mDy^yX{6t)'^, Vre > 0. 
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In order to find its global maximum, we need to compute the first derivative. 


f{w) = 6e-^o/2 


1 



\/w 


k‘l9y5t 


(1 - kyStf 1 

2ey5t jj- 


Therefore, 

f'{w) = 0 —(1 — ky5t)‘^'uP‘ + ^ywSt + kyOy = 0 . 

To solve this quadratic, we divide throughout by (St)'^ and introduce a new variable, w = w/6t. 
Then there exists a unique, positive solution tco for which the derivative is zero, namely 


Wo 


£2 


+ 


£4 

Sj/ 


+ 


k^ei 


2(1 — kydtY y 4(1 — kyStY (1 — kyStY 


From (3.31), we have that 

wq < ^-KVy ^ Wo < 2Tri>y 5t. 

Since the second root is negative, f{w) must be increasing up to wo and decreasing after this point, 
so the function attains its global maximum at wo- Hence, 

fiw) < fiwo) = < h^2TTv'^5t = rimvyiy\{5tf‘. 


Making use of the upper bound in (3.40), we derive the following inequality: 

1 + $( 2 : 0 ) - ^'(zi) < exp ^r]mvy^y\{6t)‘^'^. 


(3.41) 


Substituting back into (3.37) with (3.41), we get 


X < exp < rjm{ky6y + Uy^y)X{6t)'^ + -aXwdt >. 


Note from (3.33) that this holds when re = 0 as well. Applying (3.26) with ui = leads to 

'rf‘m‘^iyX{6t)‘^ — 2rimky5t — 2r] + 2 < 0. 


Hence, from ( |3.38[ ), 
Therefore, 


a < 2r]{m + 1 ) — 2 . 

X < exp < rjm^kyOy + Vy^y')X{5t)‘^ — Xw6t + r]{m + l)XwSt >. 


Substituting back into (|3.32|) gives the inductive step. Finally, taking m = N in (|3.30) leads to 

(3.42) 

□ 


E [0t] < exp <j ^rjXT‘^{ky9y + Uy^y) + r?ATyo 


The right-hand side is finite and independent of 6t and the conclusion follows. 

The next two results establish uniform moment bounds for both the original and the discretized 
variance and short rate processes. Note that under assumption (Al3), the moment boundedness of 
5 * and 5 *, for i £ {d,f}, extends naturally to |r*| and |f*|. 
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Proposition 3.7. The following two statements hold. 


(1) The square-root process in (3.1) has uniformly bounded moments, i.e., 

E sup yf < oo, Vp > 1. 

.te[o,r] J 


(3.43) 


(2) The FTE scheme in (3.4) for the square-root process has uniformly bounded moments, i.e., 

(3.44) 


sup E 


(5te(o,(5T) Lte[o,r] 


sup yf 


< oo, Vp > 1, V(5 t > 0. 


Proof. (1) The polynomial moments of the square-root process can be expressed in terms of the con¬ 
fluent hypergeometric function and, according to Theorem 3.1 in |24] or to they are uniformly 
bounded, i.e., 

sup E [yf] < oo, Vp >- 7^- (3.45) 

te[o,T] sj/ 

On the other hand, using Holder’s inequality, we deduce that 

sup y^ <3P~^{yo + kyeyTY+ 3P~^kPTP~^ [ yP du + sup [ y/yfdWl . (3.46) 

te[ 0 ,T] 40 i 6 [ 0 ,T] Jo 

Using the Burkholder-Davis-Gundy (BDG) inequality, we can find a constant Cp > 0 such that 

. (3.47) 


E 

sup 

/ ^dWP 

p- 

< CpE 

Y rT \p/2- 

( / Vudu] 

<\Cy+\CyTP-^¥. 

[ Vudu 


-i6[0,T] 

Jo 



\Jo J 

L 1 

.Jo 


Taking expectations in (|3.46|) and employing Fubini’s theorem, we get 
E 


sup y\ 
Lte[o,T] 


1 


< 3P-\y^ + kyOyTY + -3P-YlCy + 3P-\kl + ^.beyCy)TP sup E [yf]. (3.48) 

^ ie[o,r] 


The right-hand side is finite by (3.45), whence the conclusion. 


(2) Integrating the auxiliary process y defined in (3.3), we deduce that 


yt = yo + ky / {By - yu) du + iy / ^/yu dWP. 


(3.49) 


For any y, e > 0, there exists a constant c{p, e) > 0 such that max(0, x)^ < c{p, e)e*^^ for all x G 
In particular, this implies that yf < c(p, e)e^^‘ for all t G [0, T]. Hence, 


where 


sup E [yf] < c(p, e) sup E[Ft], 

t£[o,T] te[o,r] 


Fi = exp <j -eky / y„ du -\- ef,y / ^/yf dW^ . 

Jo Jo 


Assuming that t G [tn,tn+i] and conditioning on the cr-algebra we get 

[^t] = exp I (0.5e2^2 _ _ tn)yt„y 
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(3.50) 


(3.51) 


(3.52) 
























However, we can find e sufficiently small such that ky > 0.5e.^^. Then [Tt] < [^t^] and the 

law of iterated expectations ensures that 


E [Ft] < E [TtJ < E [Ft„_,] < ... < E [Tq] ^ sup E [Ft] = 1. 

te[o,T] 


Substituting back into (3.50), we deduce that 


sup E [yf] < c{p, e) 
te[o,r] 


(3.53) 


which is independent of 5t. On the other hand, since yt < \yt\ and using Holder’s inequality, we get 


sup yf<3P ^ {yo + kyOyTy + 3^ ^k^T^ W yPdu + 3P sup 
te[o,r] Jo ie[o,r] 


VWldwy 


(3.54) 


Taking expectations on both sides and using the BDG inequality as in (|3.47|) leads to 
E 


sup 
Lte[o,T] 


1 


< 3P-yyo + kyOyTY + -3P-yiCp + 3y-\kl + 0.5^PCp)T^’ sup E [yf], (3.55) 


te[o,r] 


□ 


for some constant Cp > 0. However, sup^gjgE [yf] < supjgjg.T] [Vt] > conclusion follows 

from (3.53). 

Proposition 3.8. Under assumption (^1), the following two statements hold. 

(1) The process gf has uniformly bounded moments, i.e., 


E 


sup (if 

(ep.TI 


< oo, Vp > 1. 


(3.56) 


(2) The process Y from (3.7) has uniformly bounded moments, i.e.. 


sup E 


(5te(o,<5T) Lte[o,r] 


sup ( 9 /)” 


< 00 , Vp > 1, V(5 t > 0. 


(3.57) 


Proof. (1) From (2.1), we know that 


9t 


^ = d + - kf f gldu-psffff a{u,Su)\l VugiduP if f \[ddWl. (3.58) 


Since gk is non-negative and using the fact that 2^/\ab\ < |o| -|- |6|, we find an upper bound 


9t — 9o kfdft -|- \psf\iff^max 

JO 


1 


1 11 / 

Vudu —( 2 kf \pgf\if(j.,r,,iax') / 9 udu-\-if I \J gL 

^ Jo Jo 

Fix t G [0,T]. Using Holder’s inequality, we get 


sup (yf)^<4P ^{gl + kfefTY + 2P ^\Psf\^Y(r?naxT^ ^ [ vPdu 
se[o,t] Jo 

+ 2P~‘^{2kf + \psf\if(TmaxYd"^~^ [ ^Up 

Jo se[o,t] 



(3.59) 
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Taking expectations on both sides and using the BDG inequality as in (3.47) leads to 


E 


sup {giY <4P + kfefTY + 2‘^P ^YfCp + 2P sup E [u^] 


- se[o,i] 


iie[o,T] 


+ (2P-Y2kf + \psf\(famaxYT^~^ + 2^P-^efCpTP-^) ^ E 
for some constant Cp > 0. Applying Gronwall’s inequality (see [E]), we get 


sup {giy 

_ se[o,M] 


du, (3.60) 


E 


sup {g{Y <Up ^{gl + kf0fTY + 2 ^p ^YfCp + 2P ^\psffYf(^?aaxT^ sup E[<]) 
V «e[o,r] ^ 


ie[0,T] 


X exp {2P-Y2kf + \psf\Cf^maxYTP + 2^P-^efCpTP} 


(3.61) 


The conclusion follows from Proposition 3.7 


(2) Integrating the auxiliary process defined in (3.6), we deduce that 


kf [ gldu-psf^f f (t{u, Su)\JVugL du + if f \fgidWl. (3.62) 


g( = gi + kfOft - 
Since g{ < \g{\, following the previous argument, we deduce that 


E 


sup {YsY <4^ ^{do +kfefTY + 2‘^P ^YCp + 2P ^\psf\^YfCr^axT^ sup E [uP] 


. se[o,i] 


ne[0,T] 


+ (2P-Y2kf + \psf\ifamaxYTP-^ + 2^P-^efCpTP-Y / E {g^Y du. (3.63) 


Since sup„g[o,T] [uS] < sup„g[o,r] Yu] and gl < sup^gjo.u] dY using Gronwall’s inequality, we get 


E 


sup {g{Y <(4P ^{gl + kf0fTY + 2^P + 2^ sup E['0P]) 

^rnT’l V ■' •' •ue[0,T] ^ 


i6[0,T] 


X exp {2P-Y2kf + MifCTmaxYT^ + 2^P-^efCpTPj . 


The conclusion follows from Proposition 3.7 


(3.64) 

□ 


Next, we establish the strong mean square convergence of the discretized variance and domestic 
short rate processes. Unlike in [32], where the focus was on the continuous-time approximation y, 
we are rather interested in the behaviour of y in the limit as the time step goes to zero. Note that 
the convergence of g'^ extends naturally to 


Proposition 3.9. The full truncation scheme in (3.5) for the square-root process converges strongly 
in the Lf sense, i.e., 

lim sup E [\yt - yY] = 0. (3.65) 

54^-0 te[o,r] 

Proof. First, fix t G [0,T]. Since \yt — yt\ < \yt — yt\-, using Cauchy’s inequality, we get 


sup \ys - ys\ < 2kyt / {yu - yu) du + 2iy sup 


sg[0,t] 


se[o,t] 


{Vyf.- Vml)dwy 


(3.66) 
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Taking expectations and then using Doob’s martingale inequality and Fubini’s theorem leads to 


E 

sup \yg 

-ys\^ 

<4klT [ E 

sup \yg - yg\‘^ 


. se[o,t] 


Jo 

.sg[0,n] 


du + AkyT"^ sup E [|yu 

iie[o,T] 


yu^ 


+ 8^yT sup E [\yu - yu\] + 8^yT sup E[|y„-y„|]. 
«e[o,r] u£[o,t] 


(3.67) 


Applying Gronwall’s inequality, we get 


E 


sup \yt - 
Lte[o,T] 







sup E 
te[o,T] 


m 


yip] +8^yT sup E [\yt 
telo,T] 


yt\] 


+ 8^yT sup E [\yt 
ie[o,r] 



(3.68) 


On the other hand, we know that 


sup E [\yt - yt\^] < 2 sup E [\yt - yt\^] + 2 sup E [\yt - yt\‘^] . 
ie[o,T] 4e[o,T] te[o,r] 


(3.69) 


Substituting into (3.69) with the upper bound in (3.68), we get 


sup E[\yt-yt\‘^]<2(l + 4klTh^^l^") sup E [\yt - yt\^] + WeyTe^'^y^" sup E[\yt-yt\] 
te[o,T] \ y J tg[o^r] te[o,r] 


+ 1642re^"^^ sup E[|yt-yi|]. 
i6[0,T] 


(3.70) 


The convergence of the three terms on the right-hand side of (3.70) is a consequence of Lemma A.3 
and Theorem 4.2 in [32], which concludes the proof. □ 


3.4 The four-dimensional system 

Even though weak convergence is important when estimating expectations of payoffs, strong con¬ 
vergence plays a crucial role in multilevel Monte Carlo methods and may be required for complex 
path-dependent derivatives. First, we prove the uniform mean square convergence of the stopped 


discretized spot FX rate process in (3.13) 


Proposition 3.10. Let Lg > Sq, Ly > vq, > g^, Lf > > k and define the stopping time 

T = inf{t > 0 : St > Lg or vt > Ly or vt > Ly or gf > or g{ > Lf or g{ < (3.71) 

Under assumptions (Ml) to (M3), the stopped process converges uniformly in Lfi, i.e., 


lim E 


(5t-s>o Lte[o,T] 


sup St/\T St/\T 


= 0 . 


(3.72) 


Proof. The absolute difference between the original and the discretized stopped processes can be 
bounded from above as follows. 


St At St At S: 


ptAr ptAr ptAr 

/ h{u){Su-Su)du + / {gi-gi)Sudu + / gt{Su-Su)du 

Jo Jo Jo 
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+ 

+ 



ctAr 


[ai- gL)Sudu + / gl{Su - Su)du 


(^a{u,Su) - a{u,Su)'^Vv^SudW^ 


+ 


+ 

tAr 


10 



Squaring both sides, taking the supremum over all t G [0, s], where 0 < s < T, and then employing 
the Cauchy-Schwarz inequality leads to the upper bound 


ns/\T f ‘1 

sup \StAr - StAr\^ < 8{Lj + L} + / (5. - Suf du + 8LlT / {gt - gifdu 

e[o,s] Jo Jo 

a{u,Su)[y^-V^jSudW^ +IQLIT / {gl - gl) du 


Note that 


+ 8 sup 
te[o,s] 

+ 8 sup 
te[o,s] 

+ 8 sup 
te[o,s] 


rtAr 


a {u, Su) y/v^ {Su - Su) dW^ 


i'sAt 


+ 16 LlT / {gl-glydu 


rtAr 


(^a{u,Su) - a{u,Su)'jVv^SudW^ 


CSAt 


ro 

{al - aiydu < / {g^^ - gUr^du < T sup Ig^^ - 

Jo iS[0,s] 


(3.73) 


(3.74) 


Hence, taking expectations in (3.73), using Fubini’s theorem, Doob’s martingale inequality and the 
Ito isometry, and upon noticing that a stopped martingale is also a martingale, we obtain 


E 


sup I 
- te[o,s] 


+ IGL^T^ sup E 


te[o,r] 


< 8L^r^ sup E 
ie[0,T] 


\dt~at\ + 8(1/^ + Lj + 4 /i^Q 3 ,)r + 32 (t^„,j,L^ 


\af-af'" 


+ ‘^‘^^liaxLlT sup E 

i6[0,T] 


\Vt - Vt\ 


E 


fsAr 


\Su-Su?du 


LJO 


+ 82LiLy E 


rsAr 


|cj(u, 5„) - a{u,Su)\‘^du 


+ IGL^T^E 


sup I^Mr-^/Arl^ 
.te[o,s] 


(3.75) 


First, we bound the last expectation on the right-hand side of (3.75) from above. From (3.58) and 
(3^, since {g^^ - g^^l < {g^^ - g^^l, we get 



Squaring both sides, taking the supremum over all t G [0, s'], where 0 < s' < s, and employing the 
Cauchy-Schwarz inequality, then taking expectations and using Doob’s martingale inequality and 
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Fubini’s theorem, we deduce that 


E 


sup <Sk)T / E [gi-gCf‘tu<T du + Sk)T / E 

te[o,s'] J JO JO 


¥u-9if 


du 


+ 



E 


^ 9u "V gL ) l'U<r 


du + 



E 




du 


+ ^Plf^'f^max^vT I E 


^ 9u V 9u ) lu<r 


du + Splfi}a‘^^^^LyT / E \gi - gi\ 

JO L 


du 


+ / E 

JO 


\Vu - Vu\ 


du + Sp^f^j-L^LfTE 


rs'/\T 


|cr(u, S„) - a[u,Su)\^du 


L JO 


However, 


(V^Z-a/^) ^ < ^^{pLr - pLrf, (3.77) 

and hence, since s' < s, 


E 


sup \gl;,^ - gl;,^\^ <{^k)T+ ‘i2i)K + 9,pl^^)al,^^L^KT) [ E sup \g{^^ - g{^ 
.ts[o,s'] J Jo LtS[ 0 ,M] 


du 




rsAr 


\cr 


{u,Su) - cr{u,Su)\‘^du 


+ 8p,Af(^maxLfT sup E 
te[0,T] 

+ {32AfT + SplfAf^l.^xLvT'^) sup E [I 5 / - 5 /lj + SkjT"^ sup E \{g{ - g{)‘^ 

te[o,T] te[o,T] 


\vt - Vt\ 


(3.78) 


Therefore, applying Gronwall’s inequality to (3.78), since s < T, we get 

xsAt 


E 


sup \g{/,^-g{^^f 
-te[o,s] 


< SplAjLvLfTE 


L JO 


|fT(u,5„) -cr(u, Su)| du 


+ ^p‘lAhliaxLfT'^ sup E flut - Ut|j + 8 A:^r 2 sup EUgl-glf 

te[o,T] J te[o,T] 

+ [ 32 ^‘jT + Splj^'ja'^ax^v'd^'^) sup E 1^/ — 5/1 

te[o,T] 


(3.79) 


where 


/3i — Sk'jT + 32^|k + Spl^Af^‘inax^vl<‘T- 
Second, assumption {A2) on the leverage function ensures that 


E 


rsAr 


|ct(u,5„) - a{u,Su)\‘^du 


< 3H^T((5t)^“ + 3H^E 


\St - Sil'^lt^rdt 


rsAr 


/ sup I S't — S't I du 

.Jo tS[0,M] 


(3.80) 


+ 3B‘^E 

where t = . Furthermore, we know that 

psAr ps 

\Su-Sufdu< / sup \St-Stfdu< / sup \StAT - StArf du. (3.81) 

Jo Jo £G[0 ,w] 


CsAt 
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Substituting back into (3.75) with (3.79) - (3.81), we deduce that 


E 

1 “ 12 
sup |5£/\7- 5*tAr| 

<P2 Te 

1 “ 12 
sup |»S£Ar 'S'iArl 

du + Ps{6t)'^°‘ + Pi sup E 

\vt - Vt\ 


.iG[0,s] 

Jo 


te[o,T] 



+ P5 sup E 

]gf-gtf 

-1- Pe sup E 

\st -9t\ 

+ Pj sup E 


ie[o,T] 


te[o,r] 


t£[0,T] 



+ /^S 


r'^ 


/ E 

\St-Siflt<r 

Jo 



dt, 


(3.82) 


where the constants /3j, for 2 < i < 8, are defined below. 


h = + 8{Lj + L) + + 96B^LIL,{1 + 

/Is = mA‘^LlL,{l + /34 = 32aL.^^T(l + ^p^efLfT^e^^^), 

/?5 = 8,LlT\ /36 = 128^^7.2(4 + , /ly = 16L2r2(l + 8k}T\P^^), 

Ps = 96B^LlL,{l + ^plf^flfT^e^^^). 

Therefore, applying Gronwall’s inequality to (|3.82[), we obtain 


E 

1 “ |2 

sup |5*iAT *5iAr| 

< P3 {St)‘^°' + Pi sup E 

\vt - Vt\ 

+ P5 sup E 

]gf-gff 


-ie[0,T] 

1 te[o,r] 


i6[0,T] 



+ Pq sup E 

\9t -9t\ 

+ Pr sup E 

\9 

{-9{f 

+ Ps 

\St-Si\Ht<r 

i6[0,T] 


te[o,r] 



Jo 



The convergence of the first term on the right-hand side as, 5t^9 is trivial, whereas the convergence 
of the next two terms is a consequence of Proposition 3.9 We now show the convergence of the Lp 
difference between gf and gf. 

Suppose that t G [tn, tn+i)- Since \g( —gl\ < 1^/ — | and using the fact that 2^/\a^ < |a| +16|, 

we can bound the difference from above as follows: 


\gt - 9tf < (kf0f6t + 0.5\psf\CfCrmaxdtvt„ + {kf + 0.5\psf\^famax)dtg{^ + dt 1 ^/ “ 

< Ak)e){5tf + p^iy^ax^Hl + { 2 kf + \psMfarnax)\6tf{g{jy4ey{w/ - w/f. 

Hence, 


sup E 
te[o,T] 


\gt-9tf < sup Eiy+Affdt sup E[g(j 


0<n<N 

+ (2fc/ + \PsMf(ymax)‘^{dtf sup E 

0<n<N L 


0<n<N 


(3.84) 


The convergence of the term on the left-hand side follows from Propositions |3.7| and 3.8 


E 


Finally, the integrand in the last term in (3.83) can be bounded from above as follows: 

pt rt 2 ^ 


li<. 

< 45tE 


= E 


/_ [gi-9L + K'^)\Sudu+ /_ a{u,Su)y/^SudW^ 


1 


t< T 


■ ft 




ft 

2‘ 

u 

{9if + {gLf + h{uf 

Sltu<Tdu 

-h 4E 

J a{u,Su)^/vuSutu<TdW^ 
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<ALl{6tf[ sup E[( 5 rf)^j+ sup E\{gl) 
I Me[o,T] Me[o,T] 


+ ^ + '^^macrL'lLydt. 


max s 


(3.85) 


The right-hand side is independent of the time t and tends to zero as (5t—)-0 from Propositions 3.7 


and 3.8 which concludes the proof. □ 

Next, we prove the uniform convergence in probability of the discretized spot FX rate process. 

Proposition 3.11. If2kf9f > and under assumptions (.4.1) to (4.3), S converges uniformly in 
probability, i.e., 

(3.86) 


lim P ( sup pt — S't I > e ) = 0 , Ve > 0. 
<5t^-o Vte[o,r] 


Proof. Fix e > 0 and note that we have the following inclusion of events, 

< w : sup |5t(a;) — 5t(a;)| > e i C i a; : sup \St{io) — St{u;)\ > e, T{io) > T 
I t&[o,T] } I te[o,T] 

U < a; : sup |5't(a;) — 5i(t<;)| > e, < T > 

I te[o,r] J 

Hence, 

< sup l^i - 5i| > e i C i sup |S'iAT - ^tArj > e i U {r < r}. 

I teo,T] J I te[o,T] J 

In terms of probabilities of events, the previous inclusion implies that 

p( sup |5'i - 5t| > e ) < P ( sup |S'iAr - *5tAT| > e ) -hP (r < T). 

Vte[o,T] / Vie[o,T] / 

The convergence in probability of the stopped process is an immediate consequence of Proposition 


(3.87) 


(3.88) 


3.10 and Markov’s inequality. Moreover, from the definition of the stopping time in (3.71), we get 


{r < r} C I sup Ui > I U i sup Dt > I U i sup gf > Ld \ U \ sup g{ > Lf\ 

I te[o,T] j I te[o,r] J I te[o,T] J I te[o,T] J 

U < sup Si > Ls > U {tk < r}, 

I iefo.Tl J 


where we define 

Tk = inf{f >0:g{< 

However, assuming that > 1, we have that 


sup Si > Ls ^ sup xt > logLs > C sup \xt\ > logL^ [>. 
te[o,r] J I te[o,T] J I te[o,T] 


(3.89) 

(3.90) 

(3.91) 


Therefore, using (3.89), (3.91) and Markov’s inequality, we find an upper bound 


<L„-iE 

sup Vt 

+ L-1E 

sup Vt 

+ L-1E 

sup gf 

+ LJ^E 

sup g{ 


.te[o,T] . 


.ie[o,T] . 


-t&[0,T] . 


Lte[o,T] J 


+ (logL,) 


sup |Xi| 
lie[0,T] 


+ P(t,, < T). 


(3.92) 
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First, recall the formula for the logarithm of the spot FX rate process, 


xt = xo + J (^g^-gl + h{u)-^a‘^{u,Su)vu^du +J a{u, Su)y/^dW^. (3.93) 

Taking the supremum over all t G [0, T] and using the fact that a < +1, then taking expectations 

and employing Doob’s martingale inequality and the Ito isometry, we get 


sup x^ 

< 1-|-|xo|-|-2/lmaa;F-|-TE 

sup gf 

-bTE 

sup g{ 

+ 4-5cr^a2-rE 

sup Vt 

.ie[0,T] 


-t&[0,T] . 


.te[o,T] . 


.te[o,T] . 


Second, we prove the almost sure positivity of the process g^ under the Feller condition, i.e., when 
2kf9f > To this end, let a = {2kf9f — and define the function F : (0,oo) i—)> M by 

F{x) = Using Ito’s formula, we have that 


E 


rTATi^ 


^{9tatJ =P{9{,)-^I a[gl) ^^^'^\kf9f - kjg^ - Psf^f(7{u, Su)y Vugl)du 

0 

TAtk 

aCfigly^'^-^^'^^dWi. (3.95) 


1 pTATK 

+ 2 ^ Jo + "du-E 


Note that 


and also that 


(^kf9f - -a{l + a)ff = -a[Akf9f - 2^ - 2a(f] = -a^ff, 


E 


10 


l«<r« du < < oo. 


so the stochastic integral on the right-hand side of (3.95) is a true martingale. Hence, 


E 


^(dTArjJ ^ ^(do) - I IE {gQ l«<r« du + akf E {gQ “ tu<r^ du 


+ a\psf\yamax^ [ vy{gl) du. 


(3.96) 


Employing Fubini’s theorem and Holder’s inequality, we deduce that 


E 


F{gyj <F{gl)-lay} I (e {gl) W,, -2a-^kAfE {gl) 


-(!+“) 


,-(!+«) 


1 +a 


-2a ^\psf\y^(Tmax SUp E 


i6[0,T] 


1 + Q 


2(l + a) 


E 


(gl) 


-(l+a) , 


1+2q 
2(l + a) 


^U<Tk 


du. (3.97) 


However, the moments of the square-root process are bounded by Proposition 3.7 Furthermore, if 
p,q G (0,1), Cp^q > 0 and c = Cp -|- Cg > 0, then for all x > 0 we have 


X — CpX^ — CqX‘^ = — (x — CX^) -|- — (x — CX*^) 

1 \p 


= CpC^-P 


C X — C i-r X 


-b CoCl-9 


_^ /_^ \1 

C 1-5X — C 1-5X 


> —CpC^-P — CqC^-1 . 
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Therefore, the integrand in (3.97) is bounded from below by a constant and thus 


E 


nskj 


<c. 


(3.98) 


for some constant C independent of k. Since g-^ has continuous paths, gl^ = k ^ and F{gl^) = 


Hence, from ( 3.98[ ) and the positivity of F, we deduce that 

¥ {t, <T) = K-^¥.\F{gl^)\r^<T 


< k""E 


< Ck 


(3.99) 


Taking the limit as )-0 in (3.88), using the upper bounds derived in (3.92), (3.94) and (3.99), and 
then employing Propositions 3.7 and 3.8, the conclusion follows from the fact that we can choose 
Lg, Ly, Ld, Lf and k arbitrarily large. □ 

Many models with stochastic volatility dynamics, including the Heston model, have the unde¬ 
sirable feature of moment instability, i.e., moments of order higher than 1 can explode in finite time 
[5]. This can cause problems in practice, for instance when computing the arbitrage-free price of 
an option with a superlinear payoff. Furthermore, establishing the existence of moments of order 
higher than 1 of the process and its approximation is an important ingredient in the convergence 
analysis (see [23]). Since the most popular FX contracts grow at most linearly in the FX rate and 
their risk-neutral valuation involves computing the expected discounted payoff, it is useful to study 
the finiteness of moments under discounting. Define R to be the discounted exchange rate process, 


Rt = 


Foexpl-y rldu-^J Su)vudu + J a{u, Su)y/v^dw{^, (3.100) 


and let R be its continuous-time approximation, 


Rt = So exp < - r^du- 


1 




a 


u,Su)vudu+ / a(u,Su)V^dWy 


(3.101) 


The next result establishes a lower bound on the explosion time of moments of order higher than 
1 of the discounted process and is an important step in proving Theorem 2.1 (Theorem 1 in |lljl. 
which, in turn, plays a key role in the calibration of the model. 

Proposition 3.12. Let a>l. Under assumptions (Ml) and (M3), if T < T*, there exists loi > a 
such that, for all u G [l,a;i), the following holds: 


sup E < oo, 
te[o,T] 

where (p{a) = a + ■\/(a — l)a, C = ^Umax and T* is as given below. 
(1) When k < (p{a)(f, 


(3.102) 


_ 


(/?(a)2C^ — k? L 2 


TT 

— -|- arctan 


(p{aYCf^ — k'^ 


(2) When k > (^(a)C; 


T* = oo. 


(3.103) 


(3.104) 
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Proof. First, we find it convenient to define a new stochastic process L by 


1 


Lt = Soexp^hmaxt- - J a^iu,Su)vudu + J a{u, Su)y/vudW., 


As Rt < Lt for all t G [0, T], it suffices to prove the finiteness of the supremum over t of 


E [Lf] =S^E 


u 


exp \ ojhmaxt + OJ J a{u, Su)^/vfL dWf^ - - J a {u,Su)vudu 


(3.105) 


(3.106) 


Second, suppose that T < T*, with T* defined in (3.103) - (3.104). If A: < (p{a)C, since (f{-) is 
increasing on [1, oo) and by a continuity argument, we can find uji > a such that, for all to G (a, wi). 


k < and T < 


TT 

— + arctan 


- 

If A: = tp{a)C, since ip{-) is strictly increasing on [l,oo) and 


k 


- A;2y_ 


(3.107) 


lim 


u)Xa+ — k?\_2 


TT 


+ arctan 




= oo, 


we can find wi > a such that (3.107) holds for all oj G (a,t<;i). Finally, if A: > ip{a)C,, by a continuity 
argument, we can find ui > a such that, for all lo G (Q!,a;i), 


k > (p{u))C. 

Third, fix a; G (a,wi) and let the functions fu},gu} ■ (l,oo) i—)■ (l,oo) be defined by 

fuj{x) = 


(3.108) 


and 


gui{x) = - l) = (a; + 1/(0; - ^ ^ (x - 1 - y^l - l/oj) 


The first function is strictly increasing, whereas the second function attains a global minimum at 
X* = 1 + ^1-1 /w. Furthermore, they are equal at x* , so 

v2 


inf max{f^{x), gu;{x)} = gu,{x*) = (w + {u - l)a;)" 

X>1 


Consider the Holder pair (p, q) satisfying q = p/{p — R) such that 

/w — 1 

p = \ + \ - and q = 1 + 


OJ 


OJ 


OJ — 1 


(3.109) 


(3.110) 


This particular Holder pair ensures an optimal lower bound T* on the explosion time. Next, define 
the quantity a = poj'^ — oj and introduce the stochastic process 


Mt=poj a{u,Su)y/^dW^ 


with quadratic variation 


{M)t=poi^ a^{u,Su)vudu. 

Jo 
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Then we can rewrite (3.106) as follows: 


E [L-] = 


'id ^(-dhmaxi 


E 


exp 




+ - a‘^{u,Su)vudu 

^ ./n 


(3.111) 


Applying Holder’s inequality with the pair (p, q) from (3.110) and taking the supremum over [0, T], 


sup E [L^] < sup E 


i6[0,r] 


te[o,T] 


exp I Mi - -{M)t 


X E 


exp \ - quj {puj - l) / Vu du 


(3.112) 


The stochastic exponential is a martingale if Novikov’s condition is satisfied, i.e., 


E 


exp<{ ^{M)t 


< E 


1 


, 2 , .2 2 


Vn, du 


< oo. 


The finiteness of the two expectations in (3.112) follows from (3.107) - (3.110) and Proposition 3.5 
The extension to the interval [l,wi) follows from Jensen’s inequality. □ 

Proposition 3.13. Let a > 1. Under assumptions (.41) and (.43), if T < T*, there exists U 2 > a 
such that, for all u G [ 1 ,^ 2 ) and 5 t < k~^, the following holds: 


sup sup E [(J?i)‘^] < 00 , 
ste(o,ST) te[o,r] 

where if (a) = a + y^(a — l)a, C = ^Omax and T* is as given below. 
(1) When k < kip{a)(^, 


(3.113) 


(2) When k > ^if(a}C, 


rj-\^ _ 


_ 


1 

(3.114) 

f{a)C - k ■ 

4fc 


(3.115) 


Proof. For convenience, define a new stochastic process L by 

rt 


Lt = 


5oexp j/imaxt - ^ y cr‘^{u,Su)vudu + j a{u, Su)V^dWf^. (3.116) 


Since Rt < Lt for all t G [0, T], it suffices to prove the finiteness of the supremum over t and dt of 


E[{Ltr] =S’^E 


exp 


^whmaxt + iuj (T{u,Su)V^dW^ -j {u, Su)vu du 


(3.117) 


Suppose that T < T*, with T* from (3.114)-(3.115). IfA:< ^(f{a)C, by a continuity argument, 
we can find a ;2 > a such that, for all w G (a, a; 2 ), 


and T < ^ 


(3.118) 
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On the other hand, if fc > by a continuity argument, we can find W 2 > a such that, for all 

LV G {a,uj2), 

1 Ak 

k > --(p{uj)C and T < . (3.119) 


ifiuifC 


Henceforth, we argue as in Proposition 3.12 and use Theorem 3.6 to deduce the finiteness of the 
supremum over t and 6t of (3.117). □ 


Higham et al. [23] proved that for a locally Lipschitz SDE, the boundedness of the pth moments 
of the exact and numerical solutions, for some p > 2, ensures the strong mean square convergence of 
the Euler-Maruyama method. The existence of moment bounds for explicit Euler approximations, 
however, remained an open problem. Recently, Hutzenthaler et al. |2S] studied SDEs with super- 
linear ly growing coefficients and proved strong and weak divergence in for all p > 1, and hence 
that the moment-bound assumption is not satisfied. To the best of our knowledge, the uniform 
boundedness of moments of order greater than 1 of discretization schemes for the Heston model 
and extensions thereof has not yet been established - this gap in the literature was also identified 
in Kloeden and Neuenkirch |3nj - and our Proposition 3.13 is the first result to address this issue. 

Since the typical payoff of a EX contract grows at most linearly in the exchange rate, it suffices 
to know the strong convergence of the discounted process in to deduce the convergence of the 
time-discretization error to zero. The following theorem can be generalized to the case relatively 
easily for all a > 1, upon noticing that the critical time T* from (3.103) - (3.104) is always greater 
than the one from (3.114) - (3.115). 

Theorem 3.14. Under assumptions (.41) to (.43), if2kf9f > 0 and T < T*, where C = icTmax 
and 

Ak , 1 




C,<2k 


+ 


c-k 


1 


the discounted process converges strongly in , i.e., 

\Rt — Rt\ 


lim sup E 
te[o,T] 


C>2k, 


= 0 . 


(3.120) 


(3.121) 


Proof. Fix e > 0 and define the event A = 

\Rt — Rt 


IR* — Rtl > e|, then 


sup E 
te[o,r] 


< sup E 
ie[o,T] 


\Rt — Rt \ i 1 a = 


-h sup E 
te[o,r] 


\Rt — Rt\^ A 


Hence, 


sup E 

i6[0,T] 


\Rt — Rt 


<e-|- sup E[Rtl^]-|- sup E[RtlA]- 
t&[o,T] t&[o,T] 


Choosing some 1 < cj < min{a;i,a; 2 } and applying Holder’s inequality to the two expectations on 
the right-hand side with the pair (p, q) = (w, returns the following upper bound: 


sup E 

i6[0,T] 


iRt — Rt 


<e + 


sup E [Rf] “ + sup E [(Rt)‘’^] ■ 
te[0,T] t£[0,T] 


sup E 
te[o,T] 


Rt — Rt\ > e 


1 - 2 - 


Note that if C ^ 2A:, then 
2 


x/C^^L2 


TT 


— -|- arctan 


> 


^/3 


> 


'C + k 


,/e-ky\ ~ ,/e-k^^c-k Q-k 
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Moreover, if A: < (" < 2k, then 




TT f k 

—h arctan — , , 

2 


2^4 /F / PX _ 4A: 


Therefore, Propositions 3.12 and 3.13 (with a = 1) ensure the boundedness of moments of order 


uj of the discounted process and its approximation. Furthermore, the convergence in probability of 
the discounted process is a simple consequence of Proposition |3.11 by taking the domestic short 
rate to be zero. Finally, taking e sufficiently small leads to the conclusion. □ 


3.5 Option valuation 

We now examine the convergence of Monte Carlo estimators for computing FX option prices when 
the dynamics of the exchange rate are governed by the Heston”2CIR++ SLV model and assumptions 
(.4,1) to (.43) are satisfied. As an aside, note that we discussed in Section]^ how other derivative 
pricing models, including popular models in equity markets, can be formulated as special cases. 
First, we consider European options. 


Theorem 3.15. Let P = E 

rT^. 


e-fo^tdt(^K_S^y 


be the arbitrage-free price of a European put 


option and P = E e -^o — St)~^ its approximation. If2kf6f > then 

lim |P — P| = 0. 

St^o 

Proof. A simple string of inequalities gives the following upper bound: 


(3.122) 


P - P < E 


11e- fo -e-fo }{K-ST)^ + e- fo - St)^ - {K - St)^}] 




- So afdt _ - /o gfdt 


+ \{K - Sry - {K - StY 


. (3.123) 


However, for any non-negative numbers x and y, \e ^ — e < \x — y\, and so we can use Fubini’s 
theorem to obtain an upper bound for the first expectation, 


(3.124) 


sup E 

1 g- fY aidu _ g- fY gidu 1 

< sup f E 

loi-dti 

du < T sup E 

Isf - 9?! 

ie[o,T] 


te[o,T] Jt 


te[o,r] 



The right-hand side tends to zero by Proposition 3.9 Define the two events A = < K} and 

A = {Sr < at}, and denote the last expectation in (3.123) by J. Then 


Therefore, 


J = E 

J < E 
< E 


(^K St) (at St) | ( Iada + ) 


|5t - StI IadA + IE {K - St) IahA' 


-hE 


{K-St) lA^nA 


{St-StIIahA +KF{AnA^)+KF{A'^nA) . 

Let 5 be an arbitrary positive number, then we have the following inclusion of events: 

4 n 4'= = ({St < a: - (5} u {a: - < Sr < a:}) n (St > k } 


(3.125) 
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C <K-6}n{ST> K}) u{K-6<St<K} 
C j St — St — (5< St < K^. 


In terms of probabilities of events, we have 

F{AnA^)<F(^\ST-ST\>6'^+F{K-S<ST<K), Vd > 0. (3.126) 

We can bound the second probability from above in a similar fashion, 

n A c I jS-r - S'rl > (5| U < 5r < K + <5} 

<P(^|5t-5t| >(^)+IP(-ftr<5r<i^ + 5), V(5>0. (3.127) 

For a suitable choice of 6, the last terms on the right-hand side of (3.126) and ( |3.127[ ) can be made 
arbitrarily small, whereas the first terms tend to zero by Proposition 3.11 Therefore, the two 
probabilities in (3.125) converge to zero as 6t —)■ 0. Finally, fix e > 0 and let 7? = {|S't — St\ > e}. 

as follows: 

\St - 5r| IahA'^b 


We can bound the expectation on the right-hand side of (3.125) as follows: 
|S’t - 5t| Iaha 


E 


< E 

< K 


|5r - S'tI 


-hE 


St — St I > e j -|- e. 


(3.128) 


Taking the limit as dt —)• 0, employing Proposition |3.H and making use of the fact that e can be 
made arbitrarily small leads to the conclusion. □ 


Theorem 3.16. Let C = E 

rT ^d. 


e- ^tdt(^ST - K)' 


and C = E 


be the arbitrage-free price of a European call 

t2 


e -^0 _ 77)^ its approximation. If 2kf9f > 0 and T < T*, with T* from 


(3.120), then 


lim \C-C\ =0. 

<5t->0 


(3.129) 


Proof. A simple string of inequalities gives the following upper bound: 

|C - C| < E [|(i?T - We" ^0^+ - {Rt - Ke--^o^^dty^ 


< KP‘ 


E 


- fo afdt _ - /g gfdt 


-7E 


|7?t — Rt\ 


(3.130) 


The two expectations on the right-hand side tend to zero as —)• 0 from (3.124) and Theorem 
3.14 , respectively. □ 


Asian options depend on the average exchange rate over a predetermined time period. Because 
the average is less volatile than the underlying rate, Asian options are usually less expensive than 
their European counterparts and are commonly used in currency and commodity markets, for 
instance, to reduce the foreign currency exposure of a corporation expecting payments in foreign 
currency. For any 0 < s <t <T, define the discount factors: 


Ds,t = e and = e 


* rfldu 


(3.131) 
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Theorem 3.17. Consider a fixed strike Asian option with arbitrage-free price 


[/ = E [e“ [ip{A{0, T) - K)]' 

[7 = E [e" ^0 T) - K)] ■ 

If2kf6f > ffj and T < T*, with T* from ( 3.120[ ), then 

lim \U -U\= 0. 

St^O 


(3.132) 


Here, Afi),T) represents the arithmetic average and f) = ±1 depending on the payoff (call or put). 
For continuous monitoring, ^( 0 ,T) = Stdt and A{0,T) = )p Stdt. 

Proof. The absolute difference can be bounded from above by 

|C/ - i7| < E [| [fi{Do,TA{0,T) - KDo^t)]^ - [^(.Do,tA(0,T) - KDo,t)Y 
Therefore, we end up with the following upper bound: 


\U-U\< Ke 


^maxT 


E 


=- So afdt _ - fi gfdt 


+ E 


\Do^tA{0, T) — Do^tA{0, T) 


(3.133) 


We deduced the convergence of the first expectation in (3.124). Using Fubini’s theorem, 
|T)o,t^(0, T) — Z)o,T^(0, r)| 


E 


< —E 

- j. 


LJO 


< sup E 

i6[0,T] 

The triangle inequality leads to the following upper bound: 

sup E \Dt^TRt - Dt;TRt\ 
te[o,T] L 


\Do^TSt — Do^TSt \ dt 


\Dt,TRt - Dt^rRtl 


< sup E 

Rt\e A 

te[o,r] 


_|_ g/lmaxT gg^p Jg 

\Rt — Rt\ 

te[o,T] 



rT .d 


du I 


(3.134) 


Since both g'^ and g'^ are non-negative processes, for any 7 greater than one we have 

- e~SA sidun < C-fA aidu _ ffi gidu\ ^ ^ [0,r]. 


Applying Holder’s inequality to the first expectation on the right-hand side of (3.134) with the pair 
(w, 7 ), where 1 < a; < wi and 7 = a;/(a; — 1 ), and using the last inequality, we find that 


sup E 

i6[0,T] 


R,\e-ft 9idu_^-fi gidn\ < g^p E [Rf]^ sup E 

ie[0,T] ie[0,T] 


- SA aidu _ - ff gidu 




(3.135) 


The convergence of the first term on the right-hand side of (3.134) is a consequence of (3.124) and 


Proposition 3.12 (with a = 1), whereas the convergence of the second term is due to Theorem 3.14 
In case of discrete monitoring or a floating strike, we follow the exact same steps. □ 
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Barrier options continue to gain popularity in many over-the-counter markets, including the FX 
market. Their popularity can be explained by two key factors. First, barrier options are useful in 
limiting the risk exposure of an investor in the FX market. Second, they offer additional flexibility 
and can match an investor’s view on the market for a lower price than a vanilla option. 

Theorem 3.18. Consider an up-and-out barrier eall with arbitrage-free price 


C/ = E 


[/ = E 


e-- K)^ t i 
^ ’ V 


™Pt6[0,T] St<B^ 

where K is the strike price and B is the barrier. If2kf0f > ffj, then 

lim \U -U\= 0. 

(5t-s.O 

Proof. Define the events ^ = { supig[o,r] ^ B^ and ^ = { sup^gjg.r] < B^, then 
\U-U\<E f| (Do,t - ^o,t) {St - K)^ t a +Do,t{ (^t - K)^ 1a-{St - K)^ Ia}] 


(3.136) 


<{B-K) 
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- /o afdt _ - fg gfdt 


+ e 


hmaxT 


E 


\{St-K)^1a-{St-K)^1a\ 


The first term tends to zero by (3.124[) and we can rewrite the second term as follows: 


E 


I {St - K)^{ Iaha^ +1 AnA) “ {St - K)^{ Iaha + ^A=nA) | 


< E 


-F E 


{St - Ky Ia^^dA +IE ISt-^tIi^ha 


(ST-i^)+l^nT 
< (B - X) + { E (A n T") + E (T" n T) } + E [|5 t - 5t| l^nA 

We can bound the last expectation from above just as in (|3.128 ) to find 


E 


|5T-5r|lAnA] <Se(|5t-5t| > e) + e, Ve > 0. 


(3.137) 


(3.138) 


Therefore, the expectation converges to zero with the time step by Proposition 3.11 Fixing <5 > 0 
and following the argument of Theorem 6.2 in |22] leads to 


An A‘^ C < sup IS"* — 5t| >(5>U<i? — (5< sup St < B 

i6[0,T] 


i6[0,T] 

In terms of probabilities of events, we have 


E(Aln^'=)<E( sup \St-St\>5)+F(B-5< sup St < 
Vie[o,T] / V te[o,T] 

We can bound the second probability in (3.137) from above in a similar fashion, 

E (a 1^ n 1) < E ( sup \St- St\>6] + f(b < sup St<B + d). 

Vte[o,T] / V te[o,T] / 

The conclusion follows from Proposition 3.11 since 6 can be arbitrarily small. 


(3.139) 


(3.140) 
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Theorem 3.19. Consider any type of barrier put option with arbitrage-free price 

fT 


U = E e" ^ {K - St) ^ 1a 
U = E\e-^o 

where the events A and A depend on the type of barrier. If2kf9f > ffj, then 

lim \U -U\= 0. 

(5t-s.O 

For instance, a down-and-in barrier is associated with the set ^ = {inftg[o,r] < B) . 
Proof. An upper bound for the absolute difference can be obtained as follows: 

I[/- C7| < E fI (Zlo.T - ^o,t) {K - St )^ 1a + ^o,t {{K - St)^ t a-{K - St)^ 1^}\ 


(3.141) 


< Ke'^ 
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- So afdt _ - fg gfdt 


+ e'‘ 


E 


\{K - St)^ 1a-{K - St)^ 1a\ 


The first term tends to zero by (|3.124[) and we can bound the second term as in (|3.137[): 


E 


{K - St)^ 1a-{K- St)^ 1 a I < E (A n + E (A" n A) | 


+ E 


{K-St)^ - {K-StY 


(3.142) 


The events A and A differ with the barrier (down-and-in, down-and-out, up-and-in, up-and-out), 
however one can show in a similar way to (3.139) and (3.140|) that 


lim E (A n A'^) = 0 and lim E (A'^ n A) = 0 
St^o St^o 


(3.143) 


for any type of barrier. Finally, the convergence of the last term on the right-hand side of (3.142) 
was derived in Theorem |3.15 which concludes the proof. 


□ 


Theorem 3.20. Consider a down-and-in/out or up-and-in barrier call option with arbitrage-free 
price 


U = E 
U = E 


e-fo ^"^*(St-K)^1a 

e-fo^^~'‘^*(ST-K)^lA 


where the events A and A depend on the type of barrier. If 2kf6f > and T < T*, with T* from 
(3.120), then 

lim |C/-^1=0. (3.144) 

St^O 

Proof. An upper bound for the absolute difference can be obtained as follows: 

\U -U\<e\\{Rt- KDo,t)^ - {Rt - KDo,t)^\ Iaha 
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{Rt - KDq^t)^ IahA'^ 


-h E 


[Rt - KDo,t)^ Ia'^ha 
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Therefore, we end up with 


VI 

1 

b 

1-1 

1 

4 


1 e“ So dUt _ g- f^ gfdt 1 

-bE 


-bE 



The convergence of the first two terms on the right-hand side is a consequence of Theorem 3.14 


where 1 < cj < min{a;i,a; 2 }, we find that: 


and (3.124), respectively. Applying Holder’s inequality with the pair (w, 7 ) to the last two terms, 

" T (A n ^ 


and 


E 


E 


Rt 1 


AnA<^ 


< E 


Rt 1 


A‘^nA 


< E 


{Rt)^^ “ p (a= n A) r. 

Using Propositions 3.12 and 3.13 (with a = 1) and the limits in ( |3.143 ) concludes the proof. □ 
Among the most developed exotic derivatives in the FX market are the double barrier options. 
Theorem 3.21. Consider a double knock-out eall option with arbitrage-free price 
U = E\e--^o^^<iCST-K)^ls- r ... 

\ ^ / vnftg[o,T] St>L, supjg[g y] St<Bjj ’ 

U = E\e-^o^^<iU St - K)^ tr - - i , 

L \ 1 J {inqg[o,T]5t>L,sup,g[o.T]>St<s| J’ 

where K is the strike and L, B are the lower and upper barriers, respeetively. If2kj9f > then 

lim|t/-U|=0. (3.145) 

St^O 

Proof. First, note that we have the following inclusion of events: 


I inf St > L, sup 5"* < n I inf St > L, sup St < BV 
*e[o,T] te[o,T] t6[o,r] te[o,T] 


te[o,r] 


i6[0,T] 


C N sup St < H} n { sup St > H} ) U f { inf 5* > L} n { inf < L| j. 


i6[0,T] 


Te[o,r] 


The rest of the argument follows closely that of Theorem 3.18 and is thus omitted. 


□ 


4 Numerical results 

In this section, we consider the pricing problem for exotic products, in particular, for structured 
notes embedding barrier features. Popular FX exotic products include the power reverse dual 
currency note (PRDC) [9], especially long-term (e.g., 30 years), and the forward accumulator [42] . 
Adding stochastic rates improves the pricing of both contracts for mid- to long-term expiries. 
Moreover, while the plain PRDC can be seen as a strip of vanilla options, adding stochastic-local 
volatility dynamics improves the pricing of some variations of the contract, like the trigger PRDC 
(a PRDC with a knock-out feature). 

Henceforth, we consider an autocallable barrier dual currency note (ABDC) that is indexed on 
the EURUSD currency pair, with USD as the domestic currency and EUR as the foreign currency, 
with quarterly coupon payments. While barrier dual currency notes are typically short-term invest¬ 
ments, we examine instead a mid-term alternative with an embedded autocallable note structure. 
Eor a nominal N (in USD), a rate (strike) K and an expiry T, the life cycle of the contract is 
described below. 
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(1) A monthly fixing schedule is defined for the exchange rate S such that: 

• If S crosses the up barrier B\jo, the contract is redeemed early and a coupon Cer is paid. 
The early redemption time is denoted by ter and can be infinite in case of no knock-out 
event. 

• If S' crosses the down barrier Bdi, a short put contract is activated at expiry. The knock-in 
time is denoted by tki- 

(2) If S is above B^i, a coupon C (in USD, in % of the nominal) is paid quarterly at coupon dates 
{ti, i = 1,... ,M}, for a total of M periods, such that the accumulated coupon (in USD) at 
expiry is 

M jjd j^d 

+^-^CERlrER<r- (4.1) 

i=l T T 

(3) At expiry, if S is below K, the down barrier has been activated and early redemption has not 
occurred, then the nominal is converted to EUR at the rate K. 


This product would suit an investor who wants to outperform the money market account by taking 
the risk of conversion of the nominal value to the foreign currency. The autocall feature lowers 
the price of the product due to early termination in case the exchange rate crosses the up-and-out 
barrier, whereas the down-and-in barrier feature provides protection against the conversion risk at 
the expense of increasing the net present value (NPV). Hence, an investor would find this product 
attractive under the current volatile market conditions if they expected a less volatile market in 
the future, which would imply a lower NPV due to the increased chance of early redemption and 
conversion. Over a longer time horizon, a very stable market would be the best outcome for the 
investor, since they would get all the coupons without converting the nominal at maturity. If 
conversion occurs at expiry, the investor receives N/K in EUR in exchange for the nominal N in 
USD. If the investor converts this amount to USD, the loss becomes N(^St/K — l). Hence, the 
profit-and-loss (PnL) of the product is 

^ jjd jjd 

^ ^^^<^ER lrER<T ~ ~ ^tki<T<ter ■ (4-2) 

i=l T' T 

Hence, this product is a yield enhancement contract with no capital guarantee. The accrued coupon 
is sometimes converted at expiry with the nominal at the rate K. Finally, the NPV of the contract 
(in % of the nominal) is 


NPV = E 


r M 

E 

. i=l 


DfCl 


ti<TER 4” '^'TeR 


Cer lrBR<r —^ lmi<T<rBR 


(4.3) 


where the expectation is taken under the risk-neutral measure. We assume the dynamics of S under 


this measure as specified in (2.1) 


Next, consider the contract parameters from Tableas well as the calibrated model parameters 
and leverage function - to the EURUSD market data from March 18, 2016 - from m- Suppose 
that barriers are monitored monthly, at the fixing dates, and coupons are paid quarterly. 

We employ the Monte Carlo simulation scheme defined in Section with 5x10^ sample paths 
and 376 time steps per year to price this contract, and the numerical results are displayed in Table 
1^ Note that we computed the percentage change in premium (NPV) with respect to the reference 
model, i.e., the Heston-2CIR-l—I- SLV model. 
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Table 3: The contract parameters of an autocallable barrier dual currency note, where the nominal N is in 
USD, the expiry T is in years, the strike K and the barriers -Buo and Bdi are in % of Bq, and the coupons 
C and Cer are in USD, in % of N. 


N 

T 

5o 

K 

Bvo 

Bbi 

C 

Cer 

100000 

5Y 

1.1271 

105% 

100 % 

95% 

2.5% 

1.5% 


Table 4: The NPV of the contract in % of N, the 95% Monte Carlo confidence interval, the early redemption 
and knock-in probabilities and the percentage change in NPV, for the autocallable barrier dual currency note 
specified in Table [funder the 4-factor hybrid SLV model (2.1), the 2-factor SLV model with deterministic 
rates and the LV model (2.8). 


Model 

NPV 

95% Cl 

ER prob. 

KI prob. 

Change 

Hybrid SLV 

1.7247 

(1.7228,1.7265) 

95.1% 

19.5% 

- 

SLV 

1.7456 

(1.7438,1.7474) 

95.1% 

19.5% 

+1.21% 

LV 

1.2079 

(1.2062,1.2097) 

95.9% 

19.5% 

-29.96% 


We infer from Table|^that the price of the contract under the 2-factor SLV model is 1.21% higher 
than under the 4-factor SLV model, whereas the early redemption and the knock-in probabilities 
are the same under the two models. This suggests that the stochastic rates have a significant impact 
on the price of the contract even for a 5-year expiry. We also infer from Table that the contract is 
highly underpriced and the early redemption probability is overestimated under the LV model (2.8). 
Furthermore, we calibrated in HD the 4-factor SLV and SV models to EURUSD market data from 
March 18, 2016, and observed an almost perfect fit to vanilla options with the former as opposed 
to a poor fit with the latter. In conclusion, we notice two things. First, pricing exotic products 
with barrier features under stochastic-local volatility dynamics is of paramount importance. It is 
a well-known fact that pure SV models underestimate the knock-out probability whereas pure LV 
models overestimate it, and the true price of the contract is believed to lie in between pure SV and 
pure LV model prices. In our case, a small difference in the probability of no knock-out, which is 
leveraged by the number of coupons detached during the life cycle of the contract, can lead to a 
big difference in the NPV, and this explains the need to use SLV models to improve the pricing 
performance. Second, adding stochastic rates to an SLV model is only relevant when pricing mid- 
to long-term structured notes, where the effect of the stochastic rates is accentuated. The impact of 
stochastic rates on the contract price is expected to increase in the presence of non-zero correlations 
between the spot FX rate and the short rates. 

A small difference in the probability of no knock-out gives rise to a big difference in the NPV 
since it is leveraged by the number of coupons detached during the lifecycle of the contract. 

We could extend the model (2.1) to multi-factor short rates when pricing exotic products where 
the rates appear explicitly in the payoff, for example, a spread option with the payoff 




Lt-K 


(4.4) 


where Lt the Libor rate at the fixing date T. In this case, the calibration algorithm described in 
Sectioncan still be applied, with a higher computational cost due to the more complex simulation 
scheme. 
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We conclude this section with an empirical convergence analysis of our Monte Carlo simulation 
scheme. Since the product can be decomposed into a linear combination of first-order exotics, the 
theoretical convergence (without a rate) of Monte Carlo estimators follows automatically from the 
analysis in Section The data in Table [^suggest a first-order convergence of the time-discretization 
error. As an aside, note that in the case of continuously monitored barriers, we could use Brownian 
bridge techniques to recover the first-order convergence. 


Table 5: The Monte Carlo estimates of the NPV of the contract specified in Table for different numbers 
of time steps (per year) and 5x 10^ sample paths (for a standard deviation of 9.29 x 10“"^), the difference 
from the previous NPV estimate and the empirical convergence order. 


Time steps 

NPV 

Difference 

Order 

12 

1.9081 

- 

- 

24 

1.8110 

0.0971 

- 

48 

1.7623 

0.0487 

0.996 

96 

1.7388 

0.0234 

1.057 

192 

1.7285 

0.0103 

1.184 


The empirical findings of this section demonstrate the importance of stochastic-local volatility 
dynamics as well as stochastic short rate dynamics for the pricing of long-dated exotic FX prod¬ 
ucts. Furthermore, we verified the convergence of our Monte Carlo simulation scheme for one such 
product. 

5 Conclusions 

Our aim was to establish the strong convergence of an Euler scheme for a hybrid stochastic-local 
volatility model. The only previous published work related to this problem that we are aware of 
is [22], which proves the convergence of an Euler discretization with a reflection fix in the context 
of Heston’s model and options with bounded payoffs. We established the strong L^-convergence 
(without a rate) of the discounted exchange rate approximation, a result which can be generalized 
to the LP case relatively easily, for all p > 1, albeit under a stronger condition on the maturity, and 
which is particularly useful when proving the convergence of Monte Carlo simulations for valuing 
options with unbounded payoffs. 

The analysis carried out in this paper can be extended to other financial derivatives, including 
digital options, forward-start options and also double-no-touch binary options, to name just a few. 
Furthermore, we may consider a multi-factor extension of the short rate model, in which case the 
convergence analysis applies with some slight modifications of the proofs. 

However, several unsettled questions remain, like the exact strong convergence rate of the full 
truncation scheme for the CIR process, or the strong convergence rate of schemes for the type of 
SDEs studied in this paper. On top of these being interesting and practically relevant questions in 
their own right, a sufficiently high order enables the use of multi-level simulation, as in m, with 
substantial efficiency improvements for the estimation of expected financial payoffs. 
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